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Abstract 



The Kerman-Klein formulation of the equations of motion for a nuclear 
shell model and its associated variational principle are reviewed briefly. It is 
then applied to the derivation of the self-consistent particle-rotor model and 
of the self-consistent cranking model, for both axially symmetric and triaxial 
nuclei. Two derivations of the particle-rotor model are given. One of these is 
of a form that lends itself to an expansion of the result in powers of the ratio 
of single-particle angular momentum to collective angular momentum, that is 
essential to reach the cranking limit. The derivation also requires a distinct, 
angular-momentum violating, step. The structure of the result implies the 
possibility of tilted-axis cranking for the axial case and full three-dimensional 
cranking for the triaxial one. The final equations remain number conserving. 
In an appendix, the Kerman-Klein method is developed in more detail, and 
the outlines of several algorithms for obtaining solutions of the associated 
non-linear formalism are suggested. 
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I. INTRODUCTION 



The aims of this paper are to study anew the foundations of the particle-rotor model 
and of the cranking model . The particle-rotor model (PRM) was introduced as an angular 
momentum-conserving phenomenological description of odd deformed nuclei. Because of its 
relative ease of application and, on the whole, quite remarkable success, it has been applied 
even up to the present (for instance [§]), with various alterations of detail, to a myriad of 
applications, over a lifetime of more than four and a half decades. Among the extensions, 
we mention in particular that to the description of triaxial nuclei [§-§1, the original model 
having been formulated for axially symmetric nuclei. 

In one of the textbooks 0, p. 109, we find, after a glowing appraisal of the success of 
the model, the following statement: "However, until now a clear-cut microscopic derivation 
has been missing." In fact, a microscopic derivation had been given earlier |Q, based on the 
Kerman-Klein (KK) method |@-r2|. The microscopic foundation of the axially symmetric 
PRM was studied more recently in [|I^ , starting from a semi-microscopic version of the KK 
approach, and compared in accuracy, for several examples of well-deformed nuclei, both with 
its more accurate progenitor and with the inherently less accurate cranking approximation. 

The cranking model was originally introduced into nuclear physics Jl^, within the frame- 
work of a prescribed single-particle model, to deal with the enigma presented by the first 



values encountered for the moments of inertia of deformed nuclei. An extended version [jT5 
the one considered in most applications until recent years, was based on the self-consistent 
mean-field theory of a deformed rotating object. This early work was designed primarily to 
provide formulas for the moment of inertia. 

The full range of applicability of the self-consistent cranking model, as well as its limi- 
tations, was realized in the so-called cranked shell model (CSM) [jl6|, that has been widely 
applied to the analysis of band-crossing and other high-spin phenomena. (For a current list 
of references, especially reviews, see fl^-) The formulations under discussion, which apply 
to axially symmetric nuclei, assume that collective rotation occurs about a principal axis 
perpendicular to the symmetry axis. Such a formulation is referred to currently as principal- 
axis cranking (PAC) as opposed to a recent generalization, called tilted-axis cranking (TAG) 
||T7|^T|?]. In the latter, even in the axial case, the system may rotate about an axis in a 
principal plane of the assumed (quadrupole) intrinsic shape, and for the triaxial case about 
an arbitrary (dynamically determined) direction with respect to the principal axes. 

A second aim of the present paper is to establish the relationship of the cranking models, 
including the recent generalized versions, to a microscopic theory. The previous literature 
on this subject is modest in extent. The standard references are ||2^j23[] , the major results 
of which are reproduced and discussed in 0. Briefly, starting from a formulation of the mi- 
croscopic theory by means of generator coordinates, the energy is evaluated approximately 
as a power series in the angular momentum by a method due to Kamlah valid for 
large deformations. When the variational method is applied to the lowest non-trivial ap- 
proximation of this procedure, it can be shown that the cranking theory is a solution of the 
resulting equations. This is summarized by stating that cranking is a solution, involving a 
semi-classical approximation, of the method of variation after projection as opposed to the 
exact procedure of variation before projection. 

To our knowledge, the only other studies of this subject are those based on the KK 
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method, a brief treatment of the case of rotation in a plane ||25| that predates the above- 
cited work and two studies that postdated them, one again on the problem of rotation in a 
plane |^ and the second a restricted study of the triaxial case |2^. (Some discussion of the 
cranking limit, also based on a variant of the KK method, can be found in p9|-) Up to now 
we have never presented a full account of the three-dimensional treatment either for axial 
or for triaxial nuclei. Since our methods contain features distinct from those found in the 
standard literature and in view of the renewed interest in generalized cranking models 
1^-^, the publication of a detailed account is perhaps justified, even at this late date. 



The foundations of the study are presented in Sec. II. We utilize a shell model Hamilto- 
nian, widely employed for medium and heavy nuclei, with two-particle interactions in which 
the latter are separated into two parts clearly distinguished as multipole and pairing forces, 
respectively. The advantage of such a model is that the (c-number) equations of motion that 
can be derived from it by the KK method are completely rigorous. It is a simplifying feature 
for the further study to recognize that these equations can be derived from a variational 
principle that we called the trace variational principle, suggested in our earliest paper §| 
and developed more fully in |^ and in |2^. This variational principle has several note- 
worthy features: (i) It is formulated for the many-body problem in the language of second 
quantization, (ii) The quantities varied are not wave functions, but rather a suitably chosen 
set of matrix elements, in our case coefficients of fractional parentage (to be discussed at the 
appropriate point of Sec. II). (iii) Rather than involving the Rayleigh-Ritz principle for one 
state at a time, the functional to be varied is the trace of energy expectation values over a 
prescribed space of states. 

[it turns out that not all aspects of our formulation are novel. Thus an incomplete version of the trace 
variational principle is to be found in one of the initial series of papers on matrix mechanics ]30|] , in which 
the variational parameters are matrix elements of the coordinates and momenta. This application to particle 
quantum mechanics was discovered and developed independently by us in several accounts of which the most 
recent is A version of the trace variational principle can, furthermore, be found is in a classic text in 
mathematical physics |32|. Here the formulation is close to standard Rayleigh-Ritz, in that the quantities 
varied are wave function. This formulation has found its way into the theory of density functionals |3^ ] and 
even been generalized to include the case where the trace is replaced by a different weighted diagonal sum 
Most recently the trace variational principle for fields has appeared in a quaternion generalization of 
quantum mechanics ||35[|.] 

The theory is elaborated in Sec. II only as far as is required for the remaining body 
of the text. Further development is presented in Appendix A, with an eye to formulating 
algorithms that can eventually be applied to the study of the self-consistent problem posed 
by the formulation in Sec. II. We turn to applications in Sec. Ill, where we derive the self- 
consistent PRM from the variational principle associated with the KK equations. (With one 
possible exception we are unaware of any recent work, other than our own, that has 



examined the foundations of the PRM.) 

The formalism presented in Sec. Ill does not lend itself naturally to a derivation of the 
self-consistent cranking theory, which should be a limit of the self-consistent PRM. In Sec. 
IV we describe an alternative derivation of the PRM, following ideas first advanced briefiy 
in [^, that does lead directly to the cranking limit. The considerations of Sees. Ill and 
IV apply to axially symmetric nuclei. Both treatments are extended to the case of triaxial 
nuclei in Sec. V. Further discussion of results and conclusions are given in Sec. VI. 
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II. EQUATIONS OF MOTION AND VARIATIONAL PRINCIPLE 



We choose a shell-model Hamiltonian in the form 

H = haO^aa + -Fa^spO^ad-ydpaS + -^G a-yl350'aO'\0'5ap ■ (2.1) 

In this standard model, the a^, aj^ are the destruction, creation operators for fermions in the 
shell-model mode a = {nljrriT) (r distinguishing neutrons from protons); F^^sp describes 
multipole forces and Ga-ySfs pairing forces. In this version, all multipolarities allowed by 
angular momentum conservation are included, though in practice we limit ourselves to the 
lowest few multipoles of each type. We shall also consistently use the summation convention, 
except when we wish to highlight some set of indices. With the help of the definitions 

Faysi3 = s^iJarriajc - mc\LM l) s p{j dmdjb - mb\LML)FacdbiL), (2.2) 
Gaysp = {jamajcmc\LML){jdmdjbmb\LML)Gacdb{L), (2.3) 
= (-1)>-™^ = ^2j, + l{j,m,j, - m,|00), (2.4) 

where {jmj'm'\LM) is a Clebsch-Gordon (CG) coefficient, the operator equations of motion 
can be obtained in the form 

[a^, H] = h'^tta + Faa'f3'l30'a'0'pO'l3' + Gaa' /Bfi'O'a'O-p'O-p, (2.5) 
, 1 2L + 1 , , 

K = ha- 2^^bab^-^, (2.6) 

[at, H] = —h'^a)^ — FpjjiaiaO^piaisa}^, — Gp^ia'a^^a'apCi^pn (2.7) 
2/ _|_ 1 

K = K + 2^--^GababiL). (2.8) 

Here, for example, a = {ja, —ma). 

To develop a dynamical scheme, we turn to the problem of obtaining equations for the 
matrix elements of Eqs. ( pl5| ) and ( |2.7| ). We designate a state of interest of an odd nucleus 
as I JMz/), where J is the total angular momentum, M is its z component, and v are the 
remaining quantum numbers necessary for unique specification of the state. Neighboring 
even nuclei are specified, correspondingly, as \IMn), referring to a heavier neighbor, and 
|/Mn), referring to a lighter neighbor. Below we shall then derive equations for the matrix 
elements, referred to as CFP (coefficients of fractional parentage) 



{JMp\aa\IMin) = VjMv{,aIMin), (2.9) 
{JMu\al\ IMin ) = Ujmu{oiI Mm). (2.10) 

We shall require the full notation when we turn to applications in the next section. For the 
formal developments of this section, we utilize a compressed notation, with 

JMu^i, IMin^n. (2.11) 

With new symbols defined and discussed below, we thus obtain the equations 
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£,V^{an) = (el - ^^)K.(«n) + Fo.^,p,p[V^1{Pn')V,{l3'n)]V,{a'n') 

+Go.^,p-p, [U*, {^'n')V, {M]U^ (« V) , (2.12) 

SMo^n) = (-e^ - E*^U,{an) - F-p-p,^,SJ:W)U,{a'n)]U0n') 

+G^^,-p,p[V,1{^n')U,{P'n)]V,{a'n'). (2.13) 

In the definitions that follow, we understand that Ei is the energy of the state \i) and 
that En and En are, correspondingly the energies of the neighboring even states, with the 
subscript standing either for the ground state, or for the lowest energy state considered, 
which for conciseness we shall continue to refer to as the ground state. We thus encounter 
the quantities 

S, = -E, + ]^{E-^ + E^), (2.14) 
e', = K-\ (2.15) 

A = ^(Eo-i^o), (2.16) 
E: = En-Eo. (2.17) 

The physical significance of the quantities defined in Eqs. (|2.14| )- (|2.17|) is evident. Si are 
the negatives of the energies of the odd nucleus relative to the ground-state energies of its 
even neighbors, e^, variously primed, are single-particle energies measured relative to the 
chemical potential A, and E* are excitation energies of the appropriate even nuclei. Finally 
in achieving the form of Eq. (|2.13| ), we have assumed that F and G are real. Given the 
Hamiltonian ( ^.1| ), Eqs. ( |^.12| ) and ( [^.13 ) are an exact set of consequences that define a 
non- linear eigenvalue problem with eigenvalue Si. The elements on the right hand sides of 
these equations define an effective Hamiltonian that will be discussed in considerable further 
detail in the course of this work. 

We display next a functional, JF, whose vanishing first variations yield the equations of 
motion, namely, 

= e'jViian)\^ - e';^\Uiian)\'' 

\ Fo.a'p'p [V; iPn') V, {(3'n)] [V* {an) V, {a'n')] 
+Go,^,p,-p p; V, i^'n)] [V: {an)UM'n')] 
-\Fp^'P'^ [U; {(3n') U, (/3'n)] [U:{an)Ui {a'n') 
-^,[|V^,(«n)|2 + f/,(«n)|2] 

-El\V,{an)\^ - El\U,{an)\^ (2.18) 
= g- S,\\V,{an)\^ + \U,{an)W (2.19) 



One verifies that the equations of motion ( |^.12| ) and ( |2.13| ) emerge, respectively, from the 
requirements 

0. (2.20) 



6Vi*{an) 6U*{an) 
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It is natural to inquire at this point if the functional T has any simple physical signif- 
icance, in particular, if it is related to a Rayleigh-Ritz principle. To answer this question, 
we evaluate the sum 

Tr(H + H) = ^[(n|H|n) + (n|H|n)]. (2.21) 

n 

The evaluation of this sum with the aim of eventually recognizing the relevant pieces of T 
requires, in addition to the standard tool of completeness, some algebraic rearrangement of 
the trace involving the lighter system, just as was necessary in the equations of motion. We 
then find that the interaction terms match exactly those in Eq. ( p.l^ ), but that the single 
particle terms do not. Instead we find 

K^ha = ha, (2.22) 

K -^K + ^^^GaUL) + \I^^FaaM ^ K- (2.23) 

We are thus tempted to replace the functional JF, as basis for the theory, by a functional 
that contains the new single-particle energies. We do not make this change because it 
destroys the simple physical significance of the Lagrange multiplier terms in Eq. (|2.18|) to 



which we next turn our attention. In practice, these extra single-particle terms are often 
ignored anyway. 

We consider then the Lagrange-multiplier terms that appear in Eq. (|2l8| ). The relevant 
question concerns the constraints that have been imposed on the variations. Since 

^ \V,{an)\^ = Y,{n\aia^\n) = (n|iV|n), (2.24) 

ia a 

where N is the number operator, we see that the excitation energies E'i enter as Lagrange 
multipliers for the conservation of nucleons in the heavier even nucleus. Similarly the term 
involving the sum over the \Ui{an)\'^ expresses (to an additive constant) the conservation of 
nucleons in the lighter system. Finally, we show that the eigenvalue £i is (no surprise here) 
a Lagrange multiplier for an appropriate normalization condition. To see this we take the 
matrix element in the state \i) of the summed anticommutator, 

^{a„, at } = fi, ^1 = ^(2j, + 1) = E (2-25) 



3a 



and thus find 



Wm{an)\' + \U.{an)\'] = l. (2.26) 

an 

Orthogonality constraints on the solutions need not be imposed, since they follow directly 
from the equations of motion. 

There is more to the story, however. We must note that Eq. ( |2.26| ) is only a sum of 
required normalization conditions. From the summed anticommutator for each level, 

5]{a„at} = fi„, (2.27) 
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we have 



man 



[\V,ian)\^ + U,ian)\' 



(2.28) 



If Eqs. ( ^.12| ) and ( ^.13[ ) described a linear eigenvalue problem, it would be impossible to 
impose the additional normalization conditions represented by Eq. ( p.28|) . For the general 
non-linear problem, there is no a priori inconsistency; the satisfaction of these conditions 
will be a part of any fully satisfactory algorithm. The form of the normalization condition 
( p.28| ) suggests, furthermore, that it may be both useful and natural to rescale the CFP, 



Vi{an) = ^2ja + lvi{an), 
Ui{an) = ^2ja + lui{an). 



(2.29) 
(2.30) 



There is considerably more to the formal theory than what has been presented thus far. 
However, we have all the tools needed for the further development in the text and thus 
relegate the additional theoretical considerations to Appendix A. 



III. DERIVATION OF PARTICLE-ROTOR MODEL AND ITS CRANKING 
LIMIT: AXIALLY SYMMETRIC CASE 

As a first illustration of the formalism presented in the previous section, we assume 
that the even (core) nuclei are in a single axially-symmetric band \IMjK), where K is the 
component of the angular momentum along the figure axis. There are at least two cases 
where it makes some physical sense to isolate a single K value, where it is the ground-state 
band with i^T = 0, or where the band has a large i^- value and we are dealing with an isomeric 
state. 

We first use rotational invariance to study the structure of the amplitudes V and U 
defined in Eqs. ( p.9|) and (|2.1CI| ), respectively. For this purpose we introduce a complete set 



of states \R) localized in the Euler angles, R = (a/?7), where a, (3 are the usual polar and 
azimuthal angles, respectively, and write 

\IMiK) = J dR\R){R\IMjK) 

= {^-^y j dR\R)D^l^{R). (3.1) 

The identification of a scalar product of many-body states with the Wigner D function is not 
a trivial statement, but is rather an essential element in the definition of the model to be 
studied. In fact, the designation |_R) for the many-body state is insufficiently detailed and 
is made more explicit by the statement 

\R) = U{R)\OK), (3.2) 

where \0K) is an axially symmetric intrinsic state spinning with angular momentum K 
about its symmetry axis, and U (R) is the unitary rotation operator in the many-body space 
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defined by the Euler angles that specify the rotation R. For such a state, we thus note the 
relation, with (a/?7) = (n7) 

U{h-^)\6K) exp{-iK-f) = U{hO)\OK). (3.3) 

The introduction of strictly localized states is, of course, an idealization that ignores the 
reality of band termination, but it is a standard approximation for well-deformed nuclei. 

[The previous discussion and that which foUows does not take into account TZ invariance, the invariance 
of the quadrupole shape under a rotation of tt about a principal axis. To include this symmetry in the 
discussion, we replace the state \IMjK) by an eigenfunction of TZ, 

\IMjK] = ^{\IMiK) + - K)}. (3.4) 

We then imitate the arguments starting on p. 8 of The task is to sort and collect the extra terms that 
appear both in the equations of motion and in formulas for one and two-particle observables. ] 

When Eq. ( |3.1| ) is substituted into the definition ( p.9|) of V, and use is made of the 
definitions to be given below, we are thereby led to the study of an amplitude such as 

{JMiy\aa\R) = {JMiy\U{R)U^\R)aaU{R)\6K) 

= J2{JMiy\U\JM'iy){JM'u\U^^aaU\bK) 

M' 

= Y: D^ZmD'i±{R)xJM'Aja^^a.K){-iy-^-^, (3.5) 

M'Ka 

where the previous manipulations have utilized the following relations and definitions (of 
which the first two are standard): 

{JM\U{R)\JM') = MiM'(^), (3.6) 
U-\R)a,^U{R) = T.^^^D(£:{R), (3.7) 

K 

{JMu\a,m\0) = i-iy+'^XJMuijm, K) (3.8) 

The phase in ( p.8|) has been introduced for algebraic convenience. 

With the help of the integral of a product of three D functions ( Eq. (4.6.2)) and the 

application of standard symmetry properties of CG coefficients (Eqs. (3.5.15) and (3.5.16) 
of the same reference) , we find for the CFP defined in Eq. (p] 



VjmAc^IMiK) = J2 



2ja + r 

\IK) 

^{-iy^+>^^XjK-.Aja^^a.K). (3.9) 

A similar analysis carried out for the amplitude U defined in ( |2.10| ) yields the result 



UjMu{aIMjK) = J2 



2ja + 1 

^^_-^y-M+ja-na+3a+m,^M^j _ M\jama) 

X{JK - Kaja'ia\IK)(j)jK-Kau{jatia,K), (3.10) 

huAjana^K) = {JMu\al_jOK). (3.11) 
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It is most succinct to base further discussion on the variational principle ( |2.18| ). We 



evaluate this expression when the core collective states are restricted to the members of a 
single band of an axial rotor, and the states of the odd nucleus are any states that can arise 
from the coupling. Returning to a full nomenclature, this calls for the identifications 



n IMiK, n IMjK , 

i-^JMu. (3.12) 

We are assuming here that there are corresponding bands in the two even nuclei that couple 
to the given odd nucleus. In the following we shall also suppress the bar and underline in 
the CFP, understanding them from context, but continue to emphasize this distinction in 
the excitation energies E*. 

We consider first the evaluation of all terms in the variational principle depending 
on |Vi(an)p, which includes some of the single-particle terms and some of the Lagrange- 
multiplier terms. Writing E*{IK) for E^, we consider in particular the combination 

- J2{Sj, + E\IK))\VjMu{aIMiK)\\ (3.13) 

For the evaluation of the CFP V , we utilize Eq. (p. 91), renormalized, however, by a factor 



of (I/vStt^) in order that the reciprocal of such factors do not appear in the answer, i. e., 
we rescale x by this factor. For the evaluation of expressions involving the CFP U , we shall 
utilize Eq. (|3.1CI|) similarly renormalized. This rescaling will be understood throughout the 
remainder of this paper. 

The rest of this section consists of a relatively detailed account of the evaluation of the 
variational sum in the "intrinsic" system. Subsequent variation will led to the self-consistent 
version of the strong coupling PRM. Toward this end, as part of the definition of an axial 
rotor, we assume that, equally for the barred and underlined quantities, 

E*{IK) -^E*{P - K^) = E*{il + il + ij - K^), (3.14) 

where we have introduced intrinsic components of the angular momentum. The arrow 
indicates the replacement of an eigenvalue by an operator. This is done by making use of 
the appropriate one of the CG coefficients, understood as a scalar product, that appear in 



Eq. (IJ), as follows 

E\IK){JK - KaJaKa\IK) = {JK - KaJaKa\E*{P - K^)\IK). (3.15) 

Furthermore, in each term of the sum / is coupled with some ja to a value of J . From the 
structure of the CFP, it follows that we may replace / by J -|- ja in Eq. ( p.l5| ) and write 
(with ja -> i), 

E*{P - K^) = E*[{J + ]f - K^] 

= E*{f'-K') + ^l + -^ + ... . (3.16) 

It is not necessary for these considerations that E{IK) have the simple form of a rotor 
spectrum, only that it be a function as indicated. The first term of Eq. ( |3.16|) may be 
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replaced by an eigenvalue E*[J{J + 1) — -ft'^)], and the second term, the Coriolis coupling 
will be evaluated below. The third and possible higher order terms will not be studied here, 
but will be included in applications. Other than the Coriolis coupling and higher-order 
terms, the contribution from Eq. ( p.l3|) and of the remaining single-particle terms takes the 
form 

Yl i.^'a-^Jy)\XJK^KA3atia,K)\^, (3.17) 

ej, = £j, - E*[J{J + 1) - K\ (3.18) 
For the further evaluation we can write in full generality 

|^ = /(J^-ir^)J., (3.19) 

where in the simplest case / is just the reciprocal of the moment of inertia. Using the matrix 
elements of the raising and lowering operators, and collecting terms, we find, in connection 
with the second term of Eq. ( p.l6|) , the Coriolis coupling term, the following sum to evaluate 

M irriaM IJ Ka^'^ 

x(-l)^>+ "°+"H(-^M,J - M\jama)\\JK - K'ja<\IK) 

X /( JK) [^(J-ir + /€j(J+X-/€„+l) 



'{J + K- Ka){J -K + Ka + l)^ {ja " t^a){3a + + 1) 
X {JK -Ka- IjaKa + l\IK) + 2{K - Ka)Ka{JK - KaJ a^^a\I K)]. (3.20) 

With the help of the standard orthonormalization conditions for CO coefficients this reduces 
to the final result for the term under study, 

--\ X*JK-KA3al^a,K)f{JK) 



X [^(J + K-Ka){J -K + Ka + l)^(ja " ^^a){3a + «« + !) 



+ ^(J -K + Ka){J+K-Ka + l)^{ja + Ka){3a - + I) 
XX JK 

+2{K - Ka)KaXJK-Kajaf^Sjafia, K)]. (3.21) 

We outline briefly the corresponding calculation of the single-particle terms associated 
with the U coefficients. Only the following change is necessary: All barred energies associated 
with the heavier of the two neighboring even nuclei are replaced by underlined energies 
associated with the lighter of the two neighbors. For the terms corresponding to Eq. p. 17] ), 
we find 
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E K+^J,)\<PjK-.A3a>^a,K)\\ (3.22) 



For the Coriolis coupling term, we find 

\ E <^*JK-.Aja^a.K)l{JK) 



X [^(J + K - Ka){J -K + Ka+ l)^ {ja ' Ka){ja + + I] 



\J -K+ Ka){J + K-Ka + l)^{ja + Ka){ja - f^a + I) 
y<4>JK-Ka+lu{jal^a - l,K) + {K - Ka)l^a4>JK-Kau{jal^a, K)]. (3.23) 

We turn to the contribution of the interaction terms. We evaluate a typical quartic 
contribution to Eq. ( p.20|) , for example 

E = \F^a'p'pV.jMAl3'IMiK)V^Mu{l3I'Mi.K) 
xVrM'u'{a'rMi,K)V},^,,AoiIMiK) 

= Faa'b'b{L){jamaja' " ma'\Lma - ma'){jb'myjb - mi,\Lmy - mi,) 

X [{2ja + l)(2ja' + l)(2j, + l)(2j,, + l)]-5 

-Kyuijb' l^b')XjK~Ktuijbl^b)XJ'K-K^,u' {ja' l^a')Xj' K -Kau' iJa l^a) 

x{IMjJmh/ — Mi\jymy){I' Mji Jmi, — Mi/\jhmh) 
x{I'MrJ'ma' - Mi,\ja'ma'){IMiJ'ma - M,|j,m,) 
x{JK — Kbijb'Hb'\IK){JK — KbjbKb\I' K) 

X{J'K - K^,Ja'Ka>\I'K){J'K - KaJaKallK). (3.24) 

To obtain this form from the corresponding term of Eq. ( |2.18D we have utilized Eqs. ( p.9| ) 
and (|2.2| ). From angular momentum conservation, we have the relations 

rrib' — Ml = rrib — Mj/, rria' — Mji = rria — Mj. (3.25) 

The next step is to perform the sums over the magnetic quantum numbers rria and rrib- 
For this purpose, we use Eq. (6.6.27) of with assists from Eqs. (3.5.14)-(3.5.17) of the 
same reference. At this point we obtain the sum (henceforth we simplify the arguments of 
the single-particle amplitudes, jaKa,K a, etc.) 

^ = li7;^,^,^(_i)i"+'^^+«.'+v+''5'+'^i>+^+^'(2L + l) 

xXJK-^b'u{b')x*K-n,Ab)Xj'K-nju'{a')xyj^_^^^,{a) 
x{JK — Kb'jb'K,b'\IK){JK — KbjbKb\I'K) 

X{J'K - Ka>]a'l^a'\l'K){j'K - Kaja^allK) 

(3.26) 



J Ja' Ja 




j jb 


jb' 










r 
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The final summations that we expect to be able to do in general are over I and /'. 
Toward this end we need to re express the 6 — j symbols that occur in Eq. ( 3.26| ) in terms 



of 6 — j symbols that depend separately on I and /'. This can be done by the application 
of Eq. (6.2.12) of We quote the value of ^^^^ results from this transformation: 

X ^ II'... 



x(2L + l)(2/" + l) 



3a' 3a L\{ Ja I J' \ j 3a' J' I' \ 

^ 3b' 3b I"]\J I" 3b' j \ J 3b I" j 

X{JK - Kb'jb'l^b'\IK){j'K - Ka3al^a\IK) 
X{JK - KbjbHb\l'K){J'K - Ka'3a'K.a'\l'K) 

^XJK-Kyu{,3b'l^b')X*jK~K,ui.3bl^b) 

X XJ' K~K^,i/' {ja' ^a')Xj' K-Kaiy' 

(ja/ta). (3.27) 

We are finally in a position to do the sums over / and /' by the application of Eq. (6.2.6) 



of 12^. This leads to our final exact result. It turns out that the evaluation of the remaining 
interaction terms parallels that just described, differing only in the coupling constants and 
the single-particle functions that occur. For economy of expression, it is this total result 
that we finally quote: 

E = E ['^^aa'b'b{L)xjK-K^,u{b')XjK-Ki,u 'K-K^,u'i0,')Xj'K~KaL''i'^) 

I"JJ'vv'ja...Ka- 

-^i^fea'b'a(^)0Ji^-«,,;.(&')0}i^-«i,z.(^)0J'i^-«„,i^'(a')0}'i^-«„i.'(a) 



X 



+G'„a'6'6Xji^-K,,i.(fe')0}/^-«,i.(^)0J'i^-K,/^'(«')x}'i^-«„^'(a)] 

(2^ + 1) f 3a' 3a L \ 

(2 J' + 1) \ 3b' 3b I" J 

X(ja - Kajb'K,b'\l"K.b' - l^a){ja' " Ka' jbK.b\l" K.b " K,a') 

x{l"Kb' - KaJK - Ky\J'K - Ka){l" Kb - Ka'JK - Kb\J' K - Ka'). (3.28) 

We now collect the results of the calculations presented in this section. The strong 
coupling limit of the functional for the axial case is given by the sum of Eqs. ( |3.17| ), ( p. 211 ) 



(iiH), (13-231) , and (|3.28| ). By varying in turn with respect to X*jK-naui.o) and we 



obtain the equations of motion 

ejuXJK-Kau{o) = ^'aXJK-nau{a) - -f{JK)^J{J + K - Ka){J - K + Ka + I) 
X V (ja - fia)ija + f^a + l)XjK-«:a-l;^(«) 



+ ,J{J-K + Ka){J+K-Ka+ 1) V(ja + ^^a){3a ' f^a + 1)XJK 

~ l^a)l^aXjK-tiai'i^) 
+ [Faa'b'biL)xj' 

K — K^V' 

{b)xJK-K^,u{a) 

+Gaa'b'b{L)xj' 
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2 J + 1 \ Jb' Jb I J 

X(ja - Kajb'K,b'\lK.b' " t^a){ja' " Ka'jbKb\If^b - ^a') 



(3.29) 



^Ju9jK-KauW 



-f_{JK)^{J + K-na){J-K + Ka + l) 



X \J {ja - f^a){ja + f^a + l)4>JK-Ka-lp{<^) 



+ y {J - K + Ka){J + K - Ka + (ja + l^a){ja - f^a + l)(j)jK-Ka+lu{o) 

+2{K 

-[Fi,a'Ua{.L)<Pj'K~^.yu'ip')<i)*j,K^^^^,{h)(pjK-K^,y{,a') 
-'rGaa'b'b{L)(j)j'K-Kyu' {b )Xj'K-k,,u' 

^2L + 1 r J,, J, L) 
2 J + 1 \ Jb' Jb I J 

X{ja - l^ajb'l^b'\Il^b' - l^a){ja' " l^a' jbl^b\I l^b - l^a') 

X{lKb' — HaJ'K — Ky\JK — Ka)(yl'^b — Ha'J'K — Kh\JK — Ka')- (3.30) 

We add the normalization conditions for the particle-rotor model that follow from 
Eq. ( g:28D . We find 



^ja + 1. 



(3.31) 



IV. ALTERNATIVE DERIVATION AND ITS CRANKING LIMIT:AXIAL CASE 

In addition to the PRM, self-consistent or otherwise, we are interested in the cranking 
theory, valid in the limit in which a single-particle angular momentum ja may be neglected 
compared to the collective angular momentum. In principle, we should be able to derive 
this limit from the form of the theory developed in Sec. II. However, the interaction terms, 
as derived, do not provide a natural pathway to the limit sought. Therefore we start anew 
in this section, but concentrate on deriving an approximate version of the PRM in which 
an expansion in {{j)/J) has been made, the main difference compared to the previous 
calculation residing in the treatment of the interaction terms. We derive an approximate 
version of the PRM and then introduce the additional approximation necessary to reach the 
cranking limit. 

For present purposes it is convenient to work in coordinate-spin-isospin space, designated 
by X. We work with amplitudes that we refer to as coordinate coefficients of fractional 
parentage (CCFP), 

VjmAxIMjK) = {JMu\^{x)\7MiK), (4.1) 
UjMuixIMjK) = {JMiy\i>Hx)\ IMjK ), (4.2) 

where iplx) is the nucleon destruction operator at the space-spin-isospin point x. In terms 
of these amplitudes, we rewrite the variational functional JF of Eq. ( p.l9| ) as 



13 



T = [e{xx') - Sj^S{x - x') - E*{IK)S{x - x')]Vjmu{x'IMiK)VJm^{xIMiK) 
[-e(xx') - £jj{x - x') - E*{IK)6{x - x')]Ujmu{xIMjK)U*jm^{x'IMjK) 

+ ^F{xx'x"x"')Vj, M'u' {x"IMiK) Vj*,M,,, {x"'I'Mr K) 

xVjM.{x'l'MpK)Vj*M,{xIMiK) 
+G{xx'x"x"')VrM'u' ix"IMjK)U*j,M,^, {x"'I'MvK) 
xUjmu{x'I'MpK)V;m,{xIMjK) 

-]^F{x"'x'x"x)Uj,M'u'{x"IMiK)VU*j,M>Ax"'I'MrK) 

xUjMu{x'I'MpK)U}j,j,{xIMjK). (4.3) 

We have set e' = e" = e and shall adhere to this simplification for the remainder of our 
presentation. To carry out the transformation to Eq. (O), we have made use of a special 
mode transformation to a basis in which e{xx') is diagonal, 

tta = ipl,{x)ip{x), (4.4) 
ea'-Paix) = e{xx')ipa{x'), (4.5) 
Fixx'x"x"') = F^^s(sVlix)if,ix')if*six")ifpix"'), (4.6) 

G{xx'x"x"') = Ga^ps(fl{x"')(f*f^{x')(f.,{x')(fa{x). (4.7) 

The major device of the present derivation is to transform from angular momentum 
eigenf unctions to eigenfunctions localized in angle space, a technique that has already been 
exploited in Sec. IV. We base the developments on expressions for the CCFP that are derived 
by the same initial transformations that led to Eqs. ( ^.9] ) and ( ^.10[ ), namely 



( VjMuixMiK) \ _ f ^j^j^iJ)* /ON ( XJM'u{Rx,K) \ / 2/ + 1 (7) . 

Xjmu{Rx,K) = {JMu\^{Rx)\6k), (4.9) 
(l)jMu{Rx,K) = {JMu\ij\Rx)\6K), (4.10) 
i:{Rx) = U-\R)^{x)U{R). (4.11) 

When we substitute Eq. (O) into the first two terms of Eq. (O), we encounter the 
restricted completeness relation 



IM 



E DZiR)DZiR')^^ = Sin - h') exp[-.K(7 - Y)]^- (4-12) 



To perform the integral over 7', we note the relation ( ^.3] ), and in the further calculation 
we utilize basic properties of the D functions, following from their definition, Eq. ( |3.6| ), the 
first one also used extensively for the interaction term to be computed below. 



M 

^mL'(0) = W/'- (4-14) 
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Finally, invoking the rotational invariance of the Hamiltonian, which for the term under 
consideration means that e{RxRx') = e{xx'), we obtain the result (recall that factors of Svr^ 
are suppressed) 

[e(xx') - £jj{x - x')]xjM,y{x')x*jMuix)- (4-15) 

We note the corresponding contributions 

[-e(xx') - SjuSix - x')](j)jMuix)(l)jM^{x'). (4.16) 

Applying Eqs. ([4.13|) and (H:-14|) to the interaction term, we find at the stage that the 
sums over IM and I'M' have been carried out, 

^ J dRdR'F{xx'x''x''')D^^]j^,,{R-^R')D^^}:^,,,{R-^R') 
X Xjm'u W)x*jm"v iR'y)Xj'M"'u' iR'x')x*j'u^v^, (Rx) . (4.17) 
Introducing the definition 

R-^R' = 7^, (4.18) 

and replacing the integral over R' by an integral TZ, we could do the integrals exactly by 
decomposing the amplitudes x i^ito irreducible tensors. We have resisted the temptation to 
do this, since a full calculation was carried out in Sec. III. It is more illuminating, as well 
as simpler, to proceed approximately by expanding TZ about the unit matrix where ever it 
appears as the argument of a x function. This brings in at each order angular momentum 
operators acting on single-particle wave-functions and therefore dimensionally is the source 
of the expansion in {{j)/J). For the interaction term the cranking limit will arise from the 
leading term of this expansion. 

With the help of the rotational invariance of the interaction, in the present instance the 
relation for example F{RxRx' Rx" Rx'") = F{xx'x"x'"), and the orthonormality relations of 
the D functions, we reach the result 

I \ F{xx'x"x"')xjmu'{x")x*jm'Ax"')Xjm'u{x)x*jmu{x): (4-19) 

^ JMM'vv' '^■J + ^ 

which is almost the cranking limit. The remaining interaction terms may be written down 
by inspection, namely 

9 /, 1 \G{.Xx'x"x"')xjMu'{.x")(p*jj^^,^,{x"')(t)jM'u{x')x*jMyi.^) 

-^F(x'Vx"x)0jMi.'(a;")0}Mv(a;"')0JJ\/'i^(a;')0jMi.(a;)]- (4-20) 

It remains for us to calculate the Coriolis coupling terms. Remarking that E*{IK) is an 
eigenvalue of D^^j^, 

E*{IK)dZ = E*{lZ - K')dZ, (4.21) 
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using completeness and integrating by parts, we reach an intermediate stage of the calcula- 
tion, in the form, 

- J dR[E*{ll^ - K2)Dl^i*,(i?)xjAf'.(i?x)]M^L,xVM'V(i?x). (4.22) 

By distribution of the derivatives, /op Jop + jop, by noting the relation 

lp{R)x{Rx) = -j{x)x{Rx), (4.23) 



and by expansion in powers of {j/J), Eq. ( [4.22| ) becomes to first order (in the intrinsic 

system) 



- J2 E*[J{J+1)-K^]\xjmAx)\ 

JMv 



2 



+ [^^M£'(^)]b^(^)XiM'.(i?x)]Mi,4i?)x}M'v(^^) + - (4.24) 

The second term, the Coriolis coupling may be evaluated further by substituting the general 
relation ( p.l9| ) and the matrix elements for intrinsic components. 



J-D^M^M' = JiJ + M')iJ-M' + l)D^Z'-i, (4.25) 



J+^a, = ^(J-MO(J + M' + l)Da,+i, (4.26) 

yielding 

E fiJ^nyiJ - M)(J + M + l)uxJM+iu{x) 

JMu ^ 

+ y{J + M)iJ-M + l)j,xJM-iuix) + MjsXJMuix)}x*jMuix)}. (4.27) 
The corresponding calculation for the other Coriolis term yields the sum 

-Y.E:[AJ+i)-K']\hMu{x)\' 

JMu 



+ E f_{JK){-^{J -M){J + M + l)j+0JAf+i.(a:) 

JMu ^ 

+ \^{J + M){J -M + l)j_(f)jM-iAx) + Mjs(j)jMu{x)}(P*jMui^)- (4.28) 



The signs in the Coriolis terms ( [4.27|) and ( [4.28|) appear to be reversed compared to those 



encountered in Eqs. ( p.21| ) and ( |3.23| ), but when due account is taken of Eq. ( |4.23| ), there is 
no inconsistency. 

We are finally ready to discuss the cranking limit. The essential observation is that once 
the expansion to leading order in {{j) / J) has been made both in the Coriolis coupling and in 
the interaction terms, the resulting approximate functional JF presents itself as a single sum 
over J. However, angular momentum is still conserved at this juncture. We lose angular 
momentum conservation by assuming that consistent with the condition {{j)/J) << 1 we 
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may identify M and K, i. e., we may neglect the angular momentum transferred to or from 
the particle, and write, furthermore, {uo defined below) 



XjkAx) pJ{u) + lx.,{x), (4.29) 
XJK±Ux) ^ ^2J{u) + IC^xUx), (4.30) 

i. e., the amplitudes differing in K from the "central value" by a unit are assumed pro- 
portional to the central amplitude (which is defined as the cranking amplitude) up to scale 
factors C=p discussed below. Similar definitions hold for the (j) amplitudes. The factor 



a/2 J + 1 is inserted for convenience, as will be evident from Eq. ( 4.33|) given below 



These assumptions suggest the following definitions of the components of the angular 
frequency (overline and underline understood) 



uj^{K) = f{JK)C^^{J tK + 1){J±K), (4.31) 
ij3 = f{JK)K. (4.32) 

The introduction of the factors dp may appear gratuitous at first sight, but it is needed, 
as will become especially evident when we treat the triaxial case, to guarantee that in 
the cranking limit the theorem that the angular velocity is proportional to the angular 
momentum is valid p8| . 

Remembering the definition (|3.18| ) and reinstating Cartesian intrinsic coordinates for the 



Coriolis coupling terms, we obtain the cranking variational expression 

[J^/{2J{uj) + 1)] = e{xx')Xu,u{x')x*coAx) - e{xx')(j)^^{x)(t)*^^{x') 
+ {^ijiXi^^{x))x*cuAx) + {^iJi(f)iv,y{x))(l)*^^{x) 

+ - F{xxx"x" ) Xujv' {x" ) X* ^, {x'" )Xu^v{x)x*^^{x) 

+G(a;x'xV")w(a:")C.'(^'")Xc..(x')C.(^) 
-^F(x"Vx"x)0^,,(a^")C.'(^"')0c..(x')0L(x) 
-eujvXu^v{x)x*^^{x) - 6^^(f)^^{x)(f)l^{x). (4.33) 

The equations of motion that follow are number-conserving, and according to the definitions 
( [4.31| ) and ( ^.32| ) allow solutions with principal axis cranking 



V. TRIAXIAL ROTOR: CORE-PARTICLE COUPLING MODEL AND 

CRANKING LIMIT 



In this section, we assume that states of interest of neighboring even nuclei can be 
described phenomenologically by a Hamiltonian 

nc = ^aj^ + ^a,,{lf,l]} + .... (5.1) 

In the calculations to be described below, we shall retain only the first term of Tic- The 
underlying model arises as follows: We assume that we can identify states of the appropriate 
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even nucleus as |/M/ncr), which we read as the nth state of angular momentum / belonging 
to a triaxial intrinsic structure a. We also define a rotated intrinsic state 

\Ra) = U{R)\6a). (5.2) 

It is part of the definition of the model that the scalar product 

{Ra\IMina) ^ Pj^^R) (5.3) 

satisfies the eigenvalue equation 

nM2n = E*iIn)Fi;,l. (5.4) 

Further useful equations satisfied by or defining the model include 

\IMina) = \IMiKa)c^K^ (5.5) 

K 

W = I] c£^*cifri (5.7) 



n 



{Ra\IMKa) = f-^D^iR)- (5-8) 



We turn to the evaluation of the terms in the variational functional JF. We shall follow 
the methods of both Sec. Ill and Sec. IV, depending on the aim of a particular fragment of 
the calculation. Starting from the representation 

VjMuialMina) = {JMv\a^\Ra)F^,]SR) 

= {JMu\ao.\Ra)D^ji,]AR)c^Kn, (5.9) 

we can derive a formula for the current version of the CFP V that is analogous to Eq. ( ^.9[ ) , 
namely 



KKa ' J<1 ' 

x{IMjJ - M\jama){JKjaf^a\IK + /€,) 

xi-ir^'^^X.JKAja^a^yfl^n. (5-10) 

■iy^"'XJKuU^<y) = {JKv\a^^\Qa). (5.11) 



The corresponding formula for the CFP U is 

UjMu{aIMina) = ^ 



KKa 



2ja + 1 

x{JKjaKa\IK + Ka)(t)jKu{jal^aCr)c^K+Kan^ (5-12) 
(t>JMu{iaK.a(y) = ( JMz/|a]^_^j0o-) . (5.13) 



With these formulas, we find the contributions of the simplest single-particle terms to 
take the form, in the shell-model or mode representation, 

^ [{^a- Sj,y)\XjKu{a)\'^ 

JKvjaKa 

-{ea + £ju\(l>jKu{a)\\ (5.14) 

We study next the term involving the Lagrange multiplier E*{In). Here it is convenient 
to carry out the calculation by a method analogous to that utilized beginning with Eq. ( [4.21D . 
With the help of the defining equation ( p. 41 ) and a subsequent integration by parts, we have 
first of all 

E*{In)VjMu{IMjna) = J dR[n,{h){JMv\a^\Ra)]F'i,';l{R). (5.15) 
With the help of the completeness relation 

E <^i(^)^iMi(^') = m - R'), (5.16) 

IMin 

we thus find for the total term 

E*{Ina)\VjMu{aIMina)\'^ = J dR[nc{ii){JMv\a^\Ra)]JMv\a^\RaY . (5.17) 

The square bracket may be reexpressed as 

HM{JMvK\Ro) = [nMD\iiiR)Dli±iR)]{JKi.\a,^jOa). (5.18) 
As far as the application of Tic in ( |5.18D is concerned, we then write 

= nc{J^) + -^J^ + ... , (5.19) 

oJi 

and work only to the order indicated explicitly. 
At the same time it is convenient to rewrite 

H,{h) = i6i(/+/- + /-/+) + h^ill + P) + (5.20) 
ai = hi + 62, a2 = hi - 62, as = ^3- (5-21) 

It is now straightforward to calculate the contributions arising from the two terms of 
Eq. ( [5.19| ). For the first term we find 

-87r2{[i6i[J(J+ 1) - K^] + h,K^]\xjKu{a)? 
1 



^h^{J - K + 2){J - K +1){J + K -1){J + K)xjk-2Mx*jkM 



+ -b2^JiJ + K + 2){J + K + l){J-K-l){J-K) 

^XJK+2u{a)x*jKM}, (5.22) 
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and for the second term, 



-S'K^\h)i^{J - K + 1){J + K){ja + Ka + l){]a " lia)XJK~lu{]al^a + 



+ ^^2\/ {J - K + 1){J + K){ja - Ka + l)(ja + 'ia)XJK~lu{ja'ia " lo") 



+ -b2J{J + K+1){J- K){ja + Ka + a ^a)Xj K+lu 

+hKKaXJK,.ijaK.aCT)]XjKuijaK.aCT). (5.23) 

Both of these terms can be identified as famihar structures. By means of this identifica- 
tion we shall have achieved both a simpler form for the particle-rotor formalism and for its 
limiting case, the cranking formalism. First consider Eq. ( ^.22|) . Note that the content of 
Eqs. ( |5.4| ) and ( |5.5| ) can be rewritten as 

'^cD^MjK = ^m]k'0^c)k'k, (5.24) 
in)KK'4rl = E*{In)4:\ (5.25) 

This eigenvalue equation was associated with even nuclei and thus with integer values of 
the angular momentum. By analytic continuation, we can define a corresponding eigenvalue 
equation for odd nuclei as follows: 

{HM))kk'C^kI = E\jT)cil (5.26) 

where J, K are now half-integral. We then see that if we introduce a new set of particle 
amplitudes XJtv by means of the equation 

XJKuijK,) = C-KrXjTv{jl^), (5.27) 

we can transform Eq. ( |5.22| ) into the form 

- E\jT)\xjru{a)\\ (5.28) 
Finally, as we did for the axial case, we can combine energy terms by means of a definition 

ej, = £j, + E*{Jr). (5.29) 

We turn our attention next to Eq. ( ^.23| ). We note first that this expression is an 
expanded version of 

- [aiJi3iXJKu{a)\x*jKM), (5.30) 

where Jj acts on the value of K and ji acts on the value of Kq. Transforming to the new 
amplitudes XJti/, expression ( ^.30| ) becomes 

- [aiJi3iXJrv{,a)]x*jrM), (5.31) 
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where now 



JiXjTy = Xjr'v{Jr'\Ji\JT), (5.32) 

{jT'\J,\Jr) = c%{JK'\J,\JK)c^j^l*. (5.33) 

For the purpose of taking the cranking hmit and comparing the forms derived in Sec. 
IV, we rewrite the results found so far and the corresponding terms involving (p amplitudes 
in coordinate space. For this we require only Eqs. ( p. 27] ), the corresponding equations 

XJKuix) = {JKu\i,ix)\6a) (5.34) 
= CkIxjtAx), (5.35) 

and the similar equations for the terms involving (p. We thus find the contributions 

-^JulXjTuix)]"^ + e{xx')Xjr,yix')x*jruix) + ai[Jiji{x)XjTuix)]x*jr,yi^) 

-Lju\<PjTuix)\'^ + e{xx')(l)jruix')(f)*j^^{x) + ai[Jiji{x)(f)jru{x)](f)*j^^ix). (5.36) 

The cranking limit of these terms may now be taken by means of the replacements that 
generahze Eqs. (|4.30| ) and ([4.29|) , 



XJr'u[X) 



V2jTTCrir')xUx), (5.37) 



and t' refers to r or any of the values coupled to r by the matrices of Jj, with Cr^r) = 1. 
This is the essential blurring of angular momentum conservation that takes us from the 
conserving particle-rotor approximation to the cranking approximation. It allows us as well 
to define the components of the angular velocity in generalization of Eq. (|4.29|) . 

= aiJ2Crir')iJr'\Ji\Jr), (5.38) 

r' 

t' = r'(r). (5.39) 

As usual, there are corresponding equations for the amplitudes 0. 
We may thus replace Eq. (|5.36| ) by its cranking limit 

(2 J + l)[-£u,u\Xi^Ax)\'^ + ^ixx')xujuix')xl^ix) + u;ijiix)xu,^ix)]xlj^ix) 
-^j.|0c^i.(a;)|^ + e{xx')(l)^^{x')(l)l^{x) + a;ji(a;)0^^(a;)]0*^(2;), (5.40) 

which is indistinguishable in form from the corresponding terms of Eq. ( |4.33| ). 

It remains for us to compute the contributions of the interaction terms. We consider 
first an exact calculation analogous to that carried out in Sec. Ill, starting from the repre- 
sentations ( p. 101 ) and ( |5.12| ) for the CFP in the triaxial case. It is straightforward to follow 
the calculations that begin with Eq. ( |3.24|) and culminate with the result (|3.28|) , as soon as 
one utilizes the orthonormality relations involving the coefficients c^)j at the first step. The 
final result is 
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I"JJ'KK'uu'ja...Ka... 

+Gaa'b'bXJK-Kyiy{b')(j)*jK'~Kbiy{b)4>J'K'~K^,iy'{a')Xj'K-Kau'{(^)] 

(2J' + 1) \> J6 /" j 

X(ja - Kajb'l^b'\l"l^b' - l^a){ja' " l^a' jbl^b\l" l^b " /^a') 

X(j"/tfe/ - K^JfST - Kh'lJ'K - Ka){rKb - Ka'JK' - Kfo] J'/ST' - (5.41) 

Superficially, the change compared to Eq. (|3.28|) is that instead of a fixed value of K, we have 
a double sum over K and K'. The same expression holds for a finite number of interacting 
K bands provided the sums are restricted correspondingly. 

Finally, we consider the calculation of the interaction term by the method of Sec. IV, 
needed to obtain the cranking limit. Here, in place of Eqs. ( [4.8| )-( ^?T0D , we utilize the forms 



UjMu{xMin) ) J > \ (t)jM'u{Rx,(r) )\ ^^^"^ ^ ^ 

X.jmu{Rx,(t) = {JMu\2p{Rx)\6a), (5.43) 
^jMu{Rx,a) = {JMu\ij\Rx)\Qa). (5.44) 



Once the full completeness relation ( ^.16|) is utilized instead of the restricted complete- 



ness relation ( ^.12| ), the calculation mimics the one carried out in Sec. IV. In terms of the 
amplitudes xjtu and (pjru, the result is 

^j^{^^(3;a;V'x''')xjr^'(a;'')x}rV'(a;''')XJr'4a;')x}ri.(a^) 

Jtt'uu' 

+G'(a;x'x"a;'")xj..'(a;")0}.'.'(^"')0Jr'.(x')x}..(x) 

-iF(x"Vx"x)0j.,,(x)0},,,,(x')0j.v(x"')0},.(x")}- (5.45) 

The cranking limit of this expression is indistinguishable from the corresponding terms 
of Eq. ( [4. 331) just as was the case for the single-particle terms ( |5.40|) . Thus the form of 
the cranking variational principle for the triaxial is indistinguishable from that for the axial 
case and need not be written again. It is understood, however, that we are dealing with 
full three-dimensional cranking, and that the single-particle wave functions have suitably 
modified symmetry. 



VI. SUMMARY AND DISCUSSION 

We have studied the microscopic foundations of the particle-rotor model and of the 
cranking model for both axial and triaxial nuclei. The microscopic model was chosen in a 
form in which the interaction is given at the outset as a sum of multipole and pairing forces. 
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We carried out the study from the point of view of the Kerman- Klein method based on 
the equations of motion for single fermion operators, and this choice of interaction has the 
advantage that the c-number equations of motion for the CFP are formally exact. These 
equations of motion and an associated variational principle, worked out in Sec. II, form the 
basis for the remaining considerations. 

[in the earliest papers on the KK approach |^,^ , a more general shell-model interaction was used in the 
derivation of c-number equations. An essential part of the derivation involved the introduction of the physical 
arguments needed to separate this interaction into multipole and pairing contributions. As a consequence 
of the limitations of this procedure, the equations of motion found from it are not exact. Nevertheless, the 
final equations are formally equivalent to those utilized in this paper. The explanation for this concordance 
is that in the approach in this paper, the "error" involved in the separation has already been built into the 
starting Hamiltonian, as a further compromise, widely accepted, in the definition of the microscopic theory.] 

As the first application, we derived in Sec. Ill a self-consistent particle-rotor model for 
axially symmetric nuclei. The derivation was carried out using basic ideas developed in 
|13| , where, starting from a semi-microscopic version of the theory, we derived the standard 



non-self-consistent version of the particle-rotor model. The present discussion complements 
the previous one in the sense that it starts from the beginning and carries the reasoning up 
to the edge of the semi-microscopic form. 

We began this work with the prejudice that, as apposed to previous treatments, a natural 
path to the cranking model involved passing through the particle-rotor model. Though we 
were ultimately able to confirm this prejudice, the version of the particle-rotor model derived 
in Sec. Ill, though a useful one for applications [jl3|, does not appear to be useful in the 
further transition to cranking. For this purpose we must be able to expand all contributions 
in powers of ((j )/</), the ratio of a characteristic single-particle angular momentum to the 
collective angular momentum. We have not discovered such an expansion for the interaction 
forms derived in this first treatment. Therefore, in Sec. IV we start anew, utilizing an 
approach already described briefly for two-dimensional rotations in an early publication 
p5| . Rather than pushing through to a formally exact result, we stop the calculation at 



the leading order of the small parameter, and thus obtain an approximate version of the 
particle-rotor model that still conserves angular momentum, but is only a step away from 
the cranking limit. This further step violates angular momentum conservation by the way 
in which an angular velocity is introduced to replace the collective angular momentum. In 
Sec. V the considerations of both previous sections are generalized to the triaxial case. 

Several special features of our treatment should be highlighted. For the axial case, as 
soon as the neighboring even nuclei are represented by bands with non-vanishing K values, 
we have tilted cranking in its simplest form. A fortiori, in the triaxial case we derive the 
possibility of full three-dimensional cranking. Within our mode of analysis, these statements 
may be taken to have the status of theorems. Another feature of our derivations of cranking 
models is that number conservation is maintained. 

Nevertheless, in the light of recent developments associated with tilted cranking |[T7|-|2T 



possible limitations on our work have to be addressed. Superficially, our results apply to one 
quasiparticle spectra of odd nuclei, whereas the current focus of interest is on at least two 
quasiparticle spectra of even nuclei, and even more on multi-quasiparticle states. In principle, 
however, these examples are covered by our considerations. Thus the two quasiparticle case 
is readily derived from the formalism developed in Appendix A. The multiple quasiparticle 
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case is covered if one replaces the reference ground states of the even nuclei by suitably chosen 
band heads of two quasiparticle bands. Details of such calculations are best addressed within 
the framework of specific applications. 



APPENDIX A: FURTHER DEVELOPMENT OF THE THEORY: ITERATIVE 

SOLUTION SCHEMES 

The theory developed in Sec. II was sufficient for the purposes of the remainder of the 
body of this work, a study of the strong-coupling limit. However, we cannot resist the 
temptation to show how this theory can be developed further, and the full architecture used 
to suggest algorithms for the solution of the non-linear problem thus defined. We point out 
before proceeding, moreover, that applications of algorithms similar to some of those to be 
described were carried out in our early work p9| , ^ . 

To begin the extended development, it is helpful to introduce a more concise represen- 
tation of the equations of motion ( p.l2| ) and ( |2.13| ) for the CFP by defining the vector 

The equations of motion can then be written 

Ei^i{an) =nian,pn')^i{pn'), (A2) 
where T-C{an, f3n') is the Hermitian matrix 



TC{an, (3n') = 

( {e'^ - El)5nn'5ap + T{an, (3n') A(an, (3n') 

\ A*(/?n', an) {-e'!, - E*)6nn'Sa0 + Uan, Pn') 



(A3) 



and the potentials are defined as 

f («n,7n') = F^^sAV*{f]n')V,{5n)] 

= F^-,Sf3nu{6n,(3n'), (A4) 
T{an,^n') = Fs^^,[U:{(3n'm6n)] 

= F-s^^^ti^il - n)n22{Sn,Pn'), (A5) 

= G,^/3^-7^l2(/3r^,5r^'). (A6) 
Here we have utilized a generalized density matrix, 7l{an, Pn'), defined as 

n{an,(3n') = ^,{an)^*{f3n') = ^ , (A7) 

satisfying the idempotent condition 

= nn. (as) 
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The structure displayed in the preceeding paragraph suggests that if we knew the po- 
tentials r, r, A, the chemical potential A, and the excitation energies E^, E*, we could view 
Eq. ( |A2| ) as a linear eigenvalue problem, with given Hamiltonian Ti, with eigenvalues Si and 
with solutions normalized according to Eq. ( p.26|) . In turn, this suggests at least the first 
elements of an iterative scheme for the solution of Eq. ( |A2D , where given the uth approxima- 
tion, 'H^'^^ to 7i, we define the {u+l)st approximation as the solution of the linear eigenvalue 
equation 

^(.+i)^(.+i) _ 

Thus, with the help of Eq. ([AT]), we can construct a (tentative, see below) (z/ + l)st approxi- 
mation, TZ^'^^^\ to the density matrix. Consequently, we can calculate a next approximation 
to the "potentials" . 

But how do we generate higher approximations to the energies, including A, of the 
neighboring even nuclei also needed for the next approximation to ?^? Actually the formalism 
provides several alternatives on how to proceed at this point. We describe first a method 
which adds the minimum number of conditions sufficient to determine the energies. This is 
to calculate, e. g.. 

En = {n\H\n) ^ En{n^''+'^) = Ei'^'\ (AlO) 
En = {n\H\n) ^ E„(7^('^+l)) = E^^+'\ (All) 

from which we can form the averages and differences needed for the calculation of 

Before we can actually proceed to the next iteration, we must take another step first 
recognized in our early work on the pairing problem [^]. The general theory requires 



that if we are to interpret the diagonal elements of H as eigenvalues, then the off-diagonal 
elements of H must vanish. The latter may be calculated at any stage of approximations 
by the same techniques used for the diagonal elements. In general, we shall find at any 
intermediate stage of the calculation that these conditions are not fully satisfied. To rectify 
this deficiency, we must therefore carry out orthogonal transformations in the spaces |?7,) 
and \n) in order to eliminate the off-diagonal elements. These are equivalent to (different) 
linear transformations of the CFP V and U. Such transformations are norm preserving. 

We have now (almost) defined a cycle of the present algorithm. To achieve a full measure 
of self-consistency requires yet an additional link in the chain of reasoning. This is because 
to this point, the vectors cannot be guaranteed to satisfy the normalization conditions 
( p.28| ). In general we shall find 

J2 ^Kan)<ili{an) = AA, (A12) 

man 

whereas full self-consistency requires A^ = 1. We can rectify the deficiency by rescaling the 
solution 

^''\an) ^vl/f)(«n), (A13) 
redefining the potentials, e. g.. 



f (an, 7^0 = F,^5^yA;A^[\/;('^)(/3n')v/'^^(5n)], (AM) 
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and similarly for other potentials. This requires us to fit in an additional set of iterations 
until the values = 1 are achieved. 

We next describe an alternative algorithm in which the equations of motion and normal- 
ization conditions are used in the same way as in the algorithm just described. However the 
transition from Hf-'^^ to Ti^'^'^^^ is done differently and requires further development of the 
formalism. The first step is to combine the equations of motion (|A2|) with their complex 
conjugate equations so as to eliminate the eigenvalues Si. We thereby obtain the following 
equations for the generalized density matrix TZ, 

= 7^(an, -fn")n{^n", (3n') - H{an, -fn")n{-fn", (3n'), (A15) 

i. e., we find the vanishing of the commutator, [R,7{\ = 0. 

Before undertaking the exposition of the algorithm, we exhibit an alternative derivation 
of Eq. ( |A15[ ) utilizing a variant of the variational principle ( |2.1(j| ), ( |2.19| ). Consider the 
functional 

V = g - e(an, f3n') [R?{^n\ an) - ^]7^(/3n', an)], (A16) 

where the "new" constraint with Lagrange multiplier matrix G is for normalization in the 
density matrix form. Since 

n{an,(3n'), (A17) 



6TZ{(3n', an) 

it follows that the variational condition applied to Eq. (|A16|) yields the equation 

n-en-n& + &n = o. (ais) 



From this condition, Eq. ([A15|) is readily derived by forming the appropriate commutator. 



However, Eq. (|A18|) contains additional information that we shall exploit below. Indeed, 
we shall consider the possibility of constructing an algorithm on the basis of Eq. ( [A 181 ) 



but first we describe one that utilizes Eq. ( |A15|) for the density matrix. Recalling the first 
algorithm described, let us imagine ourselves at the point where we have an approximation 
to the density matrix TZ that has been determined from the iterative procedure associated 
with the equations of motion (|A^), and that there remains the problem of computing the 
next approximation to the energies of the even neighbors. An alternative to the procedure 
that starts with Eq. ( |A10| ) is to note that the equation of motion in the density matrix 



form (A15) provides a sufficient set of equations to determine the excitation energies when 



supplemented by the number conservation conditions 

= J2\y^i(^n)\^-N 

ai 

= J2nu{an,an)-N (A19) 

a 

= J2\Ui{an)\'^-n + N 

ai 

= J2 ^22 (an, an)-n + N (A20) 
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Just as we did previously for the Hamiltonian, we should check that the number operator 
matrices J^aT^nicurijan') and J2a'^22{c(n,an') are diagonal and, if necessary, diagonalize 
them so as to improve the convergence of the procedure. 

So far we have suggested two possible algorithms which used the equations of motion 
attached in each case to a different method for calculating the energies of the even 
nuclei. We next propose two additional algorithms, also distinguished by one of the two 
methods of computing the energy in which, however, we replace the equations of motion 
for the \E' vectors by a condition related to the density-matrix formulation. Starting from 
Eq. ( [A18| ) we derive the pair 

TZH = 7^07^ = rnz, (A21) 

which imply that 

n = en. (A22) 

Substituting this result into Eq. ([A18|) , we obtain 

n - ^{'JZ,n} = 0. (A23) 

We describe an algorithm based on Eq. ( [A2ij| ) and the normalization condition TZ^ = QTZ. 
Suppose that in a uth approximation, we have 

(7^M))2 _Q7^H = 0, (A24) 

H^'^) - ^{n^''\n^''^} = y^''\ (A25) 

As an illustration of the method of steepest descents, we now choose, in order to improve 
the value of the density matrix 

57^(^) = -r]y^''\ (A26) 

where r] is an arbitrary small parameter that may be set to unity for 3^^^^ sufficiently small. 
We can check that this choice preserves the norm to first order, since it satisfies the required 
condition 

{n^^l (57^(")} - nsn^'^ = o, (A27) 

that follows from Eq. ( |A24[ ). To see this, note that Eq. ( |A27D implies that 

7^M57^(-)7^M = 0, (A28) 

which besides the trivial solution, is satisfied by Eq. ( A26 ), as follows from Eq. ( |A25|) . 
This as far as we can go without attacking specific models. 
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